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Abstract 

We calculate the branching ratio for B — > J/ip + anything, within the color- 
singlet approximation for J/tp production, but including perturbative QCD correc- 



tions beyond the leading logarithm approximation. Such higher order corrections 
Jv ■ are necessary, in order to obtain a result that is not strongly dependent on the 

D . renormalization scale. As in the earlier work of Bergstrom and Ernstrom, we use 

a double expansion in a s and in the small ratio of Wilson coefficients L0/L2, to 
identify the dominant terms in the decay amplitude. We complete their work by 
/\ ■ calculating all the leading order terms in this double expansion. The predicted 

branching ratio is then B(B — ► J/ifi + anything) = 0.9^ ' 3 x 10~ 3 , which is well 
below the experimental value B exp = (0.80 ± 0.08)%. This confirms the suspicion 
that non-perturbative corrections to the color-singlet approximation for J /ip pro- 
duction in B decays are important. 
PACS: 13.25.Hw, 12.38.Bx, 14.40.Gx 



1 Introduction 

The hadronic B decays into charmonium originate in the weak transition 
b — ■> qcc (with q = s, d), followed by the hadronization of the cc pair into the 
charmonium bound state. They are examples of color-suppressed hadronic 
B decays: at the weak vertex, the cc pair is not created automatically in a 
color singlet, and so it will be harder for it to hadronize into a charmonium 
state, rather than into a pair of D — D mesons. The exclusive decays, such 
as B — > K^J/ip, and the inclusive decay B — ► J/ip + anything, as well as 
analogous decays into the charmonium states ip{2S) and Xci, are well studied 
experimentally. However, color-suppressed hadronic B decays are still not 
well understood theoretically. 

We will concentrate on the inclusive decay B — > J/ip+ anything. In order 
for the cc pair to hadronize into a J/ip, we require that it forms a color- 
singlet, with spin S = 1 and no relative velocity -- this is the color-singlet 
mechanism for J/ijj production. Other cc configurations may also hadronize 
into J/ip; they appear at higher orders in an expansion in the small relative 
velocity of the J ftp constituents [|TJ [fj] . Since they contribute incoherently to 
the B — > J/ip + anything decay rate, they can be studied separately. We 
are interested in obtaining a reliable prediction for the leading color-singlet 
contribution. Comparing our prediction to the data will determine whether 
other cc configurations are indeed important. 

At first sight, the decay rate in the color-singlet mechanism, with QCD 
corrections included in the leading logarithm approximation, appears to be 
well below the experimental value. However, the result is not satisfactory, 
as it retains a strong dependence on the renormalization scale. To obtain a 
reliable prediction, it is necessary to include higher order QCD corrections. 
This was the subject of the work of Bergstrom and Ernstrom in Ref. [[|. 
These authors have shown how a clever reorganization of the higher order 
corrections can be used to identify the relevant contributions to the decay 
rate, and how this will eliminate the strong dependence on the renormaliza- 
tion scale. Here, we complete their calculation and derive the prediction for 
the B — > J/tfj + anything decay rate, in the color-singlet mechanism. 

Up to small corrections of higher order in Agc,D/m&, the inclusive decay 
B — > J/ip + anything is described by the corresponding parton decay. In the 
next section, we give the effective weak Hamiltonian for the 6-quark decay, 
with the next-to-leading order Wilson coefficients of Ref. J| . In the following 
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section, we adopt the program outlined in Ref. [[J, and obtain the leading 
terms in the b — > qJ/ip decay amplitude. Finally, we give our numerical 
results and discuss their significance. 



The effective weak Hamiltonian for the b- 
quark decay 



The terms of interest in the AB = 1 effective weak Hamiltonian 



H. 



eff 



V 2 q =s,d 



<'l 



-C (fi)O 1 + 2C 2 ^)O 8 



contain the operators 



1 = c 7m (1-7 5 )c g 7 M (l-7 5 )& 
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(2) 
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they correspond to the b — > qcc (q = s or d) transition, which occurs at tree 
level in the weak interaction. The subscripts in Oi$ designate the singlet 
(1) or octet (8) color structure of the V — A currents in those operators. 
The Wilson coefficients Co^Gu) include perturbative QCD corrections to the 
weak vertex. They depend on the renormalization scale /j,, which effectively 
separates those QCD corrections, from the QCD effects that appear in the 
matrix elements {qJ / ip\Ox^\b) . The /x dependence of the matrix elements 
should cancel that in the Wilson coefficients, so that the final result for the 
decay amplitude is independent of the renormalization scale. 

In the leading logarithm approximation (LLA), the Wilson coefficients 
are H 



with 



CoGu) 


= Lo(ji) 


= 2L+(ji) 


-L-M, 


(4) 


cm 


= L 2 {n) = - [L+(p) + L-Qi)] , 
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whereas at next-to-leading order (NLO), they are 
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The anomalous dimensions and /3-function coefficients are 
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5 is the number of active flavors. In the modified minimal sub- 

, (5) _ onn+39 



and rif 

traction (MS) renormalization scheme, Aqcd — A-irfq = 209^3 MeV, which 

corresponds to a s (M z ) = 0.118 ± 0.003 §. The quantity k± = (NDR), 

=F4 (HV) or =1=6 — 3 (DRED) is regularization scheme dependent. A similar 

scheme dependence appears in the calculation of the matrix elements of Oi t $, 

such that the final result for the decay amplitude is regularization scheme 

independent. 



3 The b — > qJ/ip decay amplitude in the a s — 
Lq/ ' L2 double expansion 



The amplitude for b — > qJ/if> (q = s or d), when the mass of the g-quark 
is neglected, can be parametrized in terms of the coefficients g\ and g 2 , that 
multiply the two possible Lorentz structures of the amplitude H 



A 



b^qj/ip 



%™ 



fj/ip 



cq 
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The decay rate is then 



r(6 -► sJ/ifj) + T(b -> dJ/rJ}) 




mjj(l -rf 



\ gi \ 2 r{\ + 2r) + \g 2 \ 2 {2 + r) - Re{g ig *)6r 
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where r = m 2 j,^/ml. The J/?/> decay constant, /j/^,, is a non-perturbative 
parameter that describes the hadronization of a color-singlet cc-pair, with no 
relative velocity and in a spin S = 1 state, into a J /if) meson. It is defined 
by 

(0\c^c\J/ip) = rnj/^ f j/ip e? JH> , (16) 

and it can be measured from the rate for the decay J /if) — > e + e~ . 

In the LLA, the decay amplitude is obtained from the tree level matrix 
element of the effective weak Hamiltonian in Eq. |l]. Only the singlet operator 
0\ contributes (as in Fig. la), and one has 



0i 



Lo{jj) , £2 = 



(17) 



This leads to a LLA branching ratio (see Fig. 3) that depends strongly on 
the renormalization scale fi (to the point where it can actually vanish, for 
/i ~ 2.5 GeV!). As pointed out by Bergstrom and Ernstrom ||, this means 
that higher order contributions to the matrix elements of the operators Oi t s, 
beyond the LLA, need to be considered. More specifically, the strong \i 



dependence in the LLA amplitude is due entirely to the Wilson coefficient 
L (/i). An expansion of this coefficient around some fixed \i ~ m\, gives 

2 
L (n) = L (m b ) + -a s (m b )L 2 (m b ) \n(fi 2 /ml) + ... (18) 

TV 

- although the /i dependence only appears at order a s , it is important, since 
Lo(m b ) I L 2 {m b ) ~ 0.34 is of similar size as a s (mb) — 0.21. From Eq. [18|, one 
can conclude that the terms in the matrix elements of C^s, which cancel 
the strong /i dependence that comes from L , are terms of order a s L 2 . Such 
higher order terms in a s appear at the 1-loop level; they are absent from the 
LLA result, where 1)8 only contribute at tree level. 

The reason higher order terms in a s become important is the existence of a 
second small quantity in the calculation — the ratio L / L 2 . Then, a s L 2 /L 
is not small compared to unity, and the expansion in a s fails. Bergstrom 
and Ernstrom advocated instead a double expansion in both a s and Lq/L 2 . 
Since these are the two small quantities in the calculation, this procedure will 
correctly identify the dominant terms in the matrix elements of Oi$. Those 
are also the terms that are needed to cancel the fi dependence that comes 
from the Wilson coefficients. Here, we will follow this program, and calculate 
all the leading order terms in the double expansion. Inexplicably, this was 
not done in Ref. 0, where the important terms of order o? s L\ in the decay 
rate were not calculated (and the estimate that was given is incorrect). 

At leading order in the a s — L /L 2 expansion for the decay amplitude, 
there is a tree level contribution from 0\ (Fig. la), proportional to Cq, and 
1-loop contributions from (9g (Figs, lb and lc), proportional to a s C 2 . From 
the results for the Wilson coefficients in the previous section, and to the order 
that we are interested, 

3n { a s (n) 

a s (M w ) - a s (fi) 
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a,Qi) v " ' ~"«° " 2 ' ' (19) 

C 2 (fi) = LM (20) 

(where we have replaced the regularization scheme dependent parameter k± 
by K = (NDR), 4 (HV) or 5 (DRED)). After a lengthy calculation, we 



obtain, for the parameters g\^ in the decay amplitude, 



9\ 



-L {fi) + a s {fi)L 2 {fi) — 
3 9n 



F 



k' — k — 6 In I — t> 
mi . 



+ 11 






a s (M w ) - a s (n) 



(21 + 4n/) 



+ 



16 + 7r + 4rln2 + 



1 

Wo 
2r 2 



2-r 



6 — 6r + r 2 , 
H m(l 



lnr 



6 — 6r 



271" - 



.92 



^(aO-MaO 



1 2r 
9vr2-r 



4-3r 



ln(l — r) 



-2- 

4 



2- 

41n2 + 
■ 3r' 



2r 
1 — r 



lnr 



Z7T- 



(21) 



(22) 



2-r v y 2-r 

Both g>i and ^2 are free of infrared divergences. The imaginary parts 
correspond to the contribution from the on-shell intermediate state with a 
cc color-octet, in the 1-loop diagrams of Fig. lb. They have been previously 
calculated in Ref. J!]]. The ultraviolet divergence in g\ has been removed using 
the same MS renormalization scheme as in the calculation of the Wilson 
coefficients M: to leading order in the double expansion, the counter-terms 
that are needed are 



H c ., 



G 
72 



^ E v A v e 



=s,d 



1e + ln(47r) 

£ 



2a s (n) 



IX 



C 2 (/x)C?i (23) 



(for d = 4 — e dimensions). The regularization of the divergence generates 
the scheme dependent term k' = -2 (NDR), 2 (HV) or 3 (DRED), in the 
expression for gi (see the Appendix, for more details on how this term is 
generated). Notice that k' is such that it cancels the scheme dependence 
from the Wilson coefficients, parametrized by k. 

As for the dependence of g\p on the renormalization scale /i, there is also 
an exact cancellation between the /i dependence of the Wilson coefficients 
and that which originates in the 1-loop matrix elements of the 0\^ operators. 
The latter is shown explicitly in the expressions for g 12 . The former can be 
obtained from the expansion of L (fi) around \x ~ m^ in Eq. [18|, and the 
analogous results for a s (fi) and L^fi): 

L 2 {fi) = L 2 {m b ) +... , a s (fjL) = a s (m b ) + . . . , (24) 
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As in Eq. [18], the terms not shown are of higher order in the double expansion 
in a s — Lq/Li2- In our final result (see Fig. 3), we have chosen to keep the 
full expressions for a s (/i) and L o ,2(/-0, rather than use the above expansions. 
This leads to a residual // dependence in the branching ratio, since higher 
order terms in Eqs. |18| and |24| have not been matched by the corresponding 
higher order contributions to the matrix elements of Ox t $. The fact that the 
residual \i dependence is small suggests that such higher order terms in the 
matrix elements are negligible; keeping only the leading order in the double 
expansion, as we do in here, is a good approximation. 

In order to obtain the branching ratio for B — > J/ip + anything, one must 
consider, in addition to b — > qJ/ift, the contribution to the inclusive decay 
from other possible parton processes. The dominant process is & — ^ qj/i/jg, 
as it contributes to the rate at order a s L\ in our double expansion, through 
the diagrams in Fig. 2. After integrating over the three-body phase space, 
the result for the decay rate is 

r(6 -► sJ/ipg) + T(b -► dJ/ifrg) = 




m 



x 



(25) 



— r(l — r)(l + 37r — 8r 2 ) — (1 — 6r)rlnr 
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with r = rrVj/y/ml. Notice that this contribution to the B — >■ J/ if) + anything 
decay rate is of lower order in the double expansion than the contribution 
from b — > qJ/i/i- However, we have checked that the three-body phase space 
suppresses the rate, so that both contributions are quantitatively of similar 
size. This allows us to neglect corrections to b — > qj/ipg and other parton 
processes such as b — > qJ/ipgg or b — > qj/ipgq'q', that appear in the B — ► 
J/ip + anything rate at order a 2 s L\ or a s L L 2 . They are of the same order 
as our calculation for the b — ► q J/ip rate, but quantitatively smaller because 
of the phase space suppression. 

4 Results 



Using the expressions for the decay rates, in Eqs. 15 and 25, together with 



the results of Eqs. |21| and ^, we obtain our prediction for the B — > J/ip + 

7 



anything branching ratio, 
B(B -► J /if) + anything) = r B £ [ r(6 ->• gJ/^) + r(6 -»■ gJ/V>#) ] • (26) 

q=s,d 

In order to eliminate the factor of m\ that appears in the decay rates, and 
so minimize the uncertainty in our result, we divide the RHS of Eq. |26| by 
the expression for the inclusive semileptonic B decay rate 0, 

B(B -> X c e-77 e ) = r B -^f_ |V cfe | 2 m* /(^) &( — ,"*) , (27) 

1927r d mfe m& 

and multiply it by the experimental result 

B(B -> X c e~z7 e ) = (10.4 ± 0.4)% . (28) 

In Eq. |27|, 

/(*) = 1 - 8^ 2 + 8z 6 - 2 8 - 24z 4 In z , (29) 

is a phase space factor, and 

h(z^) = l-^^h(z), (30) 

with § 

fc(*) ~ (vr 2 - |) (1 - zf + 2, (31) 

is a QCD correction factor. In our quantitative analysis, we take for the 
6-quark pole mass m& = 4.8 ± 0.15 GeV and mb — m c = 3.40 GeV, as in [[K]]. 
As for the J/ip decay constant, it can be determined from T(J/ip — > e + e~) = 
(5.26 ± 0.37) keV 0. Within the color-singlet approximation (now applied 
to the J/ip decay) and including perturbative QCD corrections, 



r(J/V> -► e + e ) = ^-Q 2 c a(mj/^) 2 -^- m J/v , 



47T 

3 - vc ~' " /v/ \mj/fj 



2 r 

16a s (mj/^) 



3tt 



; (32) 



with a(mjy^) = 1/133 and a s (m^) = 0.245, we find /j/^, ~ 515 MeV. Given 
the importance of the QCD corrections of order a s , it is quite possible that 



higher order corrections, not included in Eq. |32], will be significant. Because 
of this potentially large theoretical uncertainty, we choose to factor out the 
dependence on fj/^ in our final result. 



In Fig. 3, we show the predicted B — > J/ip + anything branching ratio as 
a function of the renormalization scale /i, for fi ranging between nib/2 and 
2m{,. The error band corresponds to the uncertainty in A\j^. Adding to this 
the uncertainty in m^ and that due to the residual /i dependence, we find 

£?(£-> J /V + anything) = L/^y ) x 0.9±^ x KT 3 . (33) 
(We have not included the error due to the use of the parton process to 



describe the inclusive decay of the B meson |l lfl. That error is small when 

MS 

HiS 



compared to the sensitivity of our result to the exact value of A^- and m*,). 



The central value in Eq. |33] corresponds to /x = m b = 4.8 GeV and A-^- = 209 



MeV. 

It is clear that the prediction for the B — > J/ip + anything decay rate, 
in the color-singlet approximation for J/ip production and decay, falls well 
short of the experimental value B exp = (0.80 ± 0.08)%. Although one might 
have suspected that this was so, already from the LLA prediction, that result 
could not be trusted because of the strong dependence on the renormalization 
scale (see Fig. 3). With a prediction that is much more stable in /i, it can now 
be concluded with certainty that new, non-perturbative, contributions to the 
description of the J ftp bound-state are indeed important [|TJ . A calculation 
of the B —>■ J/ip + anything decay rate which attempts to include such 
contributions is given in Ref. 0. Our result can be used [0 to improve 
the color-singlet part of the B — » J/ip + anything decay rate that appears in 
there. 
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of Energy (T.T.) and a grant from the National Science Foundation (J. M.S.). 

Appendix 

The origin of the scheme dependent term k', in the expression for gi (see 
Eq. ^), deserves a brief comment. We summarize the more detailed discus- 
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sion of Ref. M (notice that our result for k' can be recovered from the more 
general case presented there). 

The k' term in g^ originates from the regularization of the divergence in 
the un-renormalized 1-loop b — > qJ/ip amplitude. That divergence is of the 
form 



.4 



d/iv 



-% V <K°**&i 



X 



- - 7s + ln(47r) 



(qJ/m\b) 



(34) 



with 



n 



xc[Yl a Y(l - 7s) - 7 M (1 - 75)7 tt 7 v ]c 



(35) 



In d = 4 dimensions, Oi t s form a complete basis under QCD corrections, and 
Q is proportional to 0\. Q = —24 0\. In d = 4 — e dimensions, however, the 
Cis basis must be extended to include evanescent operators E, and 



Q = a 1 {e)0 1 + E . 



(36) 



The evanescent operators E do not exist in 4 dimensions, and they do not 
contribute to the physical amplitude. On the other hand, the terms of order 
e in a\(e) will contribute to the finite part of the b — > qJ/ip amplitude, when 



inserted in Eq. |35|. In order to determine these terms, we must completely 
fix the regularization scheme by giving the form of the evanescent operators 
E (and the same form for these operators must be used in the calculation 
of the Wilson coefficients). One way to do this is to define the evanescent 
operators by the condition || 



E 



a/3,^S^- pi*- 8a 







(37) 



where the two pairs of fermion fields were removed from E, and each one was 
replaced by T = 7 M (1 + 75). Applying this operation to both sides of Eq. [36], 
we obtain 



oifej 



^00,15*^/3^ 8a (O 



1 afrfSTpfTSa) 



■\2{2-K'e + 0{e 2 )) . 



(38) 



This is the origin of the scheme dependent k' term in the expression for g\ 
of Eq. P 
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(b) 



(c) 

Figure 1: Feynman graphs contributing to B 
order in the a s — L Q /L 2 double expansion. 




qJ/i/j (q — s, d), to lowest 
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Figure 2: Feynman graphs contributing to B — > qJ/ipg (q = s, d), to lowest 
order in the a s — L / L 2 double expansion. 
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Figure 3: The B — > J/^> + anything branching ratio as a function of the 
renormalization scale \x. The theoretical curves, in LLA and in the a s — Lo/L2 
double expansion, correspond to fj/^ = 515 MeV, m^ = 4.8 GeV and A 



20911 M eV. 
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